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Abstract 

The possibility to observe vacuum electron-positron pair creation due to a powerful optical laser 
pulse is discussed. We employ a quantum kinetic formulation of the problem with a source term 
describing the vacuum pair production in a homogeneous electric field with arbitrary time dependence 
(dynamical Schwinger mechanism). For a periodic field weak in comparison with the critical value 
Ecr ~ m^/|e| , the electron-positron plasma density changes approximately periodically with twice 
the field frequency. Under these conditions, the mean value < n > for the density per period in the 
volume is a more appropriate characteristic quantity than the residual density rir which is taken 
over an integer number of field periods and calculated using the imaginary time method. The value 
< n > is proportional to the squared field intensity and does not depend on the frequency. We show 
that in terms of the parameter < n > an optical laser can be more effective than a X-ray one. We 
expect that it is possible to observe the vacuum creation effect not only by means of the planned 
X-ray free electron lasers but already at present-day optical lasers. 

1 INTRODUCTION 

QED is considered as the most advanced physical theory, many of its predictions have been proven 
experimentally with highest available precision. Nevertheless, some questions are discussed till now, e.g., 
the vacuum pair creation effect by a classical electric field P . A complete theoretical description of this 
effect has been obtained El EI E] , but there is still no experimental proof. The main problem is the 
high value of the critical electric field strength, necessary to be reached for the pair creation, namely 
Ecr = 1.3 X lO^^V/ cm for electron-positron pairs. According to the Schwinger formula, the pair creation 
rate in a constant electric field is 

r^exp -nTT— (1) 



d^xdt 47r3 ^ n2 ^ V E 

and therefore exponentially suppressed when E <^ Ecr ■ Fortunately, the situation changes qualitatively 
if the field acts a finite time only [21 El El IE] ■ In this case, the Schwinger formula as well as its analog for 
a monochromatic field (Brezin-Itzykson formula 0) become inapplicable. 

There are a few examples for physical situations where the Schwinger effect can be observed, e.g. 
relativistic heavy ion collisions jTUI, neutron stars I12| and focussed laser pulses ^21- The structure 
of a real laser field is too complicated for the analysis, because the Schwinger effect is non-perturbative 
and it requires the exact solution of the dynamical equations. That is why the approximation of the 
homogeneous electric field is used in most cases. According to different estimates [HI El El the effect 
of vacuum creation can not be observed with the presently achieved level of laser power, see also 

The recent development of laser technology, in particular the invention of the chirped pulse amplifica- 
tion method, has resulted in a huge increase of the light intensity in a laser focal spot |17l [T5| . The most 
advanced lasers produce pulses with intensities reaching 10^^ W/cm ^ and the pulse duration decreasing 
down to few oscillation periods. As the construction of X-ray free electron lasers XFELs ^Hl is now 
planned, the possibility of the experimental proof of the Schwinger effect attracts attention again. The 
non-stationary effects become important under conditions of short pulses. We use in our work the kinetic 
equation approach, which allows us to consider the dynamics of the creation process taking into account 
the initial conditions ,6j. Compared to the other treatments, the approach within the framework of a 
transport equation contains some new dynamical aspects, such as longitudinal momentum dependence of 
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the distribution functions and non-markovian character of the time evohition. It takes into account the 
effects of the field switching and statistics, as well [201 • This approach has been applied already to the 
periodical field case '7' with near-critical values of the field strength and X-ray frequencies. In particular, 
it was shown that there is an " accumulation" effect when the intensity of the field is about half critical: 
the average density of pairs grows steadily with the increase of the field period numbers. 

In the present work, we consider the other region of field parameters really achievable nowadays in 
the optical lasers: E ^ Ecr and v <^m , where v is the laser field frequency. We suggest to use in the 
criterion for the creation efficiency the mean value < n > for the density per period in the volume A'^ 
is a more appropriate characteristic quantity than the residual density rir ■ The latter is taken over an 
integer number of field periods and calculated using the imaginary time method. The main result is that 
optical lasers can generate a greater density < n > than X-ray ones. 

The work is organized as follows. Section 2 contains the statement of the problem and the necessary 
information about the kinetic equation which is used for the description of vacuum pair creation. We 
solve this equation numerically for the conditions of the SLAC experiment |21| and study some features of 
pair production dynamics. We compare here our results obtained on the non-perturbative basis with the 
predictions of another approach and show that optical lasers can be effective generators of electron- 
positron pairs during the action of a laser pulse. In Section 3, the low density approximation is considered. 
It allows to get some analytical results and to make simple estimates. Finally, in Section 4, we discuss 
some possibilities of direct experimental verification of pair production by high power optical lasers. 



2 THE KINETIC EQUATION APPROACH 

In the kinetic approach 6 , the basic quantity is the distribution function of electrons in the momentum 
representation /(p, t) . The kinetic equation for this function is derived from the Dirac equation in an 
external time-dependent field by the canonical Bogoliubov transformation method 2 : , or by the help of the 
oscillator representation j72j . This procedure is exact but valid only for the simplest field configurations, 
e.g., the homogeneous time dependent electric field with the fixed direction 

E(t) = (0,0,S(t)), E{t)^-A{t), (2) 

where the vector potential is given in the Hamiltonian gauge — (0, 0, 0, A{t)) and the overdot denotes 
the time derivative. Such a field is not appropriate for a quantitative description of the laser pulse, but 
can probably be used as qualitative model to estimate results. The corresponding kinetic equation in the 
collisionless limit has the form 6 



i A(p, t) J dt' A(p, t') [1 - 2/(p, t')] cos 0(p, t',t), (3) 



to 



A(P,0 - ^^W^^^, (4) 



rf/(p,t) 
dt 

where 



w(p,t) = ^m-2+pl + [p^^-eA{t)f, (5) 
t 

e{p,t',t) = 2 J dticj{p,h). (6) 
t' 

Eq. can be transformed to a system of ordinary differential equations, which is convenient for a 
numerical analysis 

/ = iA.r, 
Vl - A(l - 2/) - 2ww2, 

V2 — 2ujvi, (7) 
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Figure 1: Time dependence of the density n{t) in the volume A'^ in the weak periodic field with the 
parameters E/Ecr — 4.6 • 10^^ and v/m = 4.29 • 10~^ , corresponding to the SLAC experiment PT] 
(solid line) and in the near-critical field case of an X-ray laser with E/ Ecr = 0.24 and v/m = 0.0226 
(dashed hne). 



where vi, V2 are real auxiliary functions. The system (UJ is integrated via the Runge-Kutta method with 
the initial conditions /(p, to) = ''^i(P7 ^o) — t'2(Pj ^o) = . The momentum dependence of the distribution 
function is defined by means of a discretization of the momentum space in a 2-dimensional grid, where 
the system ^ is solved in each of its nodes. The concrete grid parameters depend on the field strength, 
where typical values are Ap « 0.05 m (grid step) and Pmax ~ (5 — 10) m (grid boundary). The particle 
number density can be found after that as a moment of the distribution function 

"(*) = '/(0^(P'*)- 
Let us consider a harmonic time dependence of the field 

Em 

E{t) = E,n smvt, A = ^ cos vt. (9) 

The time dependence of density for the field ^ with the parameters Em/ Ecr = 4.6 • 10^^ and v/m ^ 
4.29 • 10~^ , corresponding to SLAC experiments (21] is shown in Fig. 1 in comparison to the planned 
X-ray laser JHI with E/Ecr — 0.24 and v/m = 0.0226. The pair density oscillates with twice the 
frequency of the laser field. The density value , which is evaluated in the imaginary time method 
[SJ I15| . corresponds to an integer number of field periods, = n{t = 27r/z/) , and it is negligible in 
comparison with the density value nm corresponding to the electric field maximum, Um — n{t — n /2v) . 
The mean density per period < n > is of the same order as Um ■ For the conditions of the SLAC 
experiment the ratio of < n > /ur ~ 3 • 10^^ . As a consequence, in spite of the fact that the residual 
density for the X-ray laser exceeds the one for the optical laser by a large factor, the situation is different 
regarding the mean density: the optical laser can produce more pairs per volume lambda cubed than the 
X-ray one. 



3 LOW DENSITY APPROXIMATION 

The low density approximation / ^ 1 can be used in the weak field limit E <C Ecr ■ In that case it is 
possible to obtain analytic estimates for the residual ( rir ) and the maximal ( Um ) densities the Eq. (O . 
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Figure 2: The transverse (left panel) and the longitudinal (right panel) momentum dependence of the 
distribution function at t = -k j^v (maximum of the field) for the two values of the field strength: = 
4.6 • Vd~^Ecr (solid lines) and = 3 • 10~^£'cr (dashed lines) with optical frequency v — 4.26 • 10^^ m . 



The particle density in this approximation is 



and can be transformed to [JHl 



to 
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Let us assume additionally that the condition 
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is satisfied, where 7 is the adiabaticity parameter [S]. This relation can be treated as the condition 
for quasi-classical charge transport in an external field on the time scale ~ 1/^^ , if only the pairs are 
created with vanishing momenta. The latter condition was often used in relation to the longitudinal 
momentum |24j but the real momentum distribution of the electron-positron pairs has a width of the 
order of the inverse mass for both transverse and longitudinal momenta, see Fig. 12 The momentum 
distribution shape varies essentially at the moments of time corresponding to the field minima: a complex 
quasi-periodic structure with a mean period of about the inverse laser frequency is formed. The mean 
period value of such a structure decreases proportionally to the number of field periods. 
By means of the inequality (|12|l Eq. Ullfl is reduced to 

2 

1 



n{t) 



2(27r)a 



f2 



j dti eE{ti) exp{2iujti 



(13) 



The time integral is calculated analytically for the field jnj and to — with the result 



nit) = 



2(27r)3 



eE„ 



4^2 



i^^(l + cos^ lyt) + 4cj^ sin^ i^t 



2i/[i/ cos j/t cos 2cji + 2(jj sin i^t sin 2ijji] >. (14) 
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According to Eq. 114|l the residual pair density after N periods is rir — n{2TTN/v) and mean pair 
density per period < n> are estimated as 



(2 



<n> = 



4(27r) 



2 ^ eE 



- 4cj2 



/47riVw\ 
l_cos(^-— j 



Now we omit the fastly oscillating term in ()15(l and suppose additionally that 



Then we obtain the simple estimate 



(n) 



(15) 
(16) 

(17) 
(18) 



The mean density of electron-positron pairs is defined in this case only by the field amplitude and does 
not depend on the frequency within a wide range of parameters (|12|l . After the integer period number 
(when the electric field vanishes) the overwhelming part of pairs is absorbed and the residual density, 
which is estimated within the usual approach jl5| . is negligible in comparison with the mean one used 
above. For the Terawatt Nd-glass laser with the wavelength 527 nm and the field strength Em — 6 • 10^° 
V/cm |21) we have m/v w 2 • 10^ , so the mean density exceeds the residual one by more than 10 orders 
of magnitude. According to Fig. 1, there are « 10^ pairs in a volume of wavelength cubed on the average 
for one period of the laser field. The same pair density can be achieved under the conditions of an X-ray 
laser [J]. Let us notice, that the formula from Ref. 9 for the pair creation probability with the condition 

m 



w 



{eEf 



eE 
2mLU 



gives only a negligible creation probability ~ 10 in this case, 
formula l(TI?)) is not applicable for field pulses of finite duration. 



(19) 

This is not surprising because the 



4 SUMMARY 

The simplest laser field model (jHl predicts the existence of a dense electron-positron plasma during a laser 
pulse duration, which is absorbed almost completely after switching off the field. The mean density is 
defined by the field strength and does not depend on frequency. The plasma density reaches 10^^ cm ~^ 
within the range of the really achieved fields 10^" — 10^^ V/cm. The usual recipe for an experimental 
proof of the Schwinger effect suggests that an increase of the residual plasma density can be achieved 
by an increase of both the frequency and the power of the laser radiation, e.g., via XFEL facilities jl9| . 
It is possible that the simpler way is to try to probe the plasma with some external field in addition to 
the generating one. SLAC experiments j21j can serve as an example of this type, where a high energy 
( 46 GeV) electron beam crosses the focus of the 527 nm Terawatt pulsed Nd-glass laser and 106 ± 14 
positrons above background have been observed. The observed positrons are interpreted as arising from 
a two-step process: laser photons are backscattered to GeV energies by the electron beam and a collision 
between the high-energy photon and several laser photons produces an electron-positron pair. 
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